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^v I In higher-curvature inflation models (R + a n R n ), we study a parametric preheating of a scalar 

^^ , field x coupled non-minimally to a spacetime curvature R (£Rx 2 )- In the case of _R 2 -inflation model, 

^\ ' efficient preheating becomes possible for rather small values of £, i.e. |£|< several. Although the 

maximal fluctuation \/(v 2 ) ~ 2 x 10 17 GeV for £ ~ —4 is almost the same as the chaotic inflation 

V \' v ' max 3 

model with a non-minimally coupled \ field, the growth rate of the fluctuation becomes much larger 
and efficient preheating i s real ized. We also investigate preheating for R 4 model and find that the 
maximal fluctuation is \/ (y 2 ) fa 8 x 10 16 GeV for £ ^ —35. 
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I. INTRODUCTION 
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The inflation is one of the most promising models for the early stage of the Universe in modern cosmology B . It 



not only gives a natural explanation for the horizon, flatness, and monopole problem but also provides us density 
perturbations as seeds for a large scale structure in the Universe B. So far, there are several models for inflation, 
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i.e. a new inflation proposed by Albrecht and Steinhardt Q, and Linde 0], a chaotic inflation by Linde [pi, and 



higher-curvature inflation by Starobinsky B. The first two models require a scalar field which is called inflaton to 
drive inflation. The chaotic inflation scenario would be natural in the sense that initial conditions for inflation are 
not restricted so much as compared with a new inflation. On the other hand, the third model is based on adding 
the higher-curvature terms to the Einstein-Hilbert action. This model is remarkable in the sense that we have an 
inflationary solution without an inflaton scalar field. That is, inflation can be realized by purely gravitational coupling. 
Several authors Jij studied some constraints of coupling with higher-curvature and investigated its various aspects 
such as the cosmic no-hair theorem. 
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An inflationary period will end when kinetic energy of an inflaton dominates over its potential energy. At this stage, 
the inflaton field begins to oscillate coherently around the minimum of its potential. The energy of the inflaton field 
will eventually be transferred to radiation, and this process is called reheating. The original version of the reheating 
scenario was considered in Ref. m just after the new inflation model is proposed. This old scenario is based on 
perturbation theory, adding a phenomenological decay term to the inflaton equation. Since the reheating temperature 
is constrained to be relatively small such as T r < 10 9 GeV and a production of GUT scale gauge boson is kinematically 
impossible, the GUT scale baryogenesis does not work well in this scenario. 

Recently, it has been pointed out that an inflaton decay will begin in a much more explosive process, called preheating 
H Q before the perturbative decay. At this stage, the fluctuation of scalar particles can grow quasi-exponentially 
by parametric resonance. The importance of this non-perturbative stage was recognized in ref. [|E)| in the context of 
the new inflation scenario. As for the chaotic inflation scenario, Kofman, Linde, and Starobinsky first investigated 
the preheating structure with both massive inflaton and a self- interacting potential pCJ ], In the case with a self- 
interacting potential A<fi /4 of a massive inflaton, the inflaton particles are produced by parametric resonance due to 
the self coupling of the inflaton field. In this model, the evolution of inflaton quanta was investigated OJl2| by using 
a closed time path formalism [05- 17 and also studied by fully non-linear lattice simulations [18]. If we adopt the value 



of A ~ 1CP 12 which is constrained by density perturbations, it was found that the fluctuation (<5</> 2 ) ~ ICU 7 Afp L will be 
produced in the preheating stage. In the case of massive inflaton without a self-coupling, however, we do not expect a 
parametric resonance. Hence one may need another scalar field \ coupled to the inflaton field, where coupling makes 
a resonant production of the x-particle possible. Especially when the coupling constant g gets large, a parametric 
resonance will turn on in the broad resonance regime. Consequently, although the expansion of the Universe reduces 
the amplitude of inflaton field adiabatically, the fluctuation of the \ field increases quasi-exponentially. As for the 
%-particle production, several numerical works have been done by making use of Hartree mean field approximations 
|fl9| and fully non- linear calculations (gFJ,|lJ. Analytic investigations including the backreaction and rescattering 
effects are performed in ref. |22|,|23| . Although the structure of resonance shows a complicated feature which is called 
stochastic, the analytic estimations based on the Mathieu equation are found to be reliable to analyze the preheating 
dynamics. 

In the higher-curvature inflation model, a reheating process in the old scenario was considered by Suen and Anderson 
p4[ . In the case of the R 2 inflation model, an effective scalar field that appears by conformal transformation can be 



approximately written as a massive scalar field in the reheating phase. Therefore, the preheating stage is absent, since 
the inflaton fluctuation cannot be amplified by a parametric resonance. In the massive inflaton model, however, a new 
scenario was recently proposed, in which the non-minimally coupled \ field is enhanced by a parametric resonance 
p5| . Especially when the coupling £ is negative, the fluctuation of the x-particle grows rapidly due to negative 
coupling instability. The final fluctuation (x 2 ) ~3x 10~ 4 A/p L in the case of £ w —4 is larger than the minimally 
coupled massive inflaton case with jkIx I0~ 3 p6| . However, in this model, rather large values of |£| are required 
to have a comparable growth rate with previous existing preheating models. Then a natural question arises. When 
one considers the non-minimally coupled x field in the context of the higher-curvature inflation model, how is the x 
fluctuation amplified? If it turns out to be larger as compared with the previous cases, it would affect the baryogenesis 
in GUT scale [£7| , a non-thermal phase transition p8| and topological defect formation |2jJ . We show in this paper 
that in the higher-curvature inflation model, especially in the R 2 inflation model, significant particle production is 
expected only due to a non-minimal coupling. As compared with the previous model of chaotic inflation with the 
non-minimally coupled scalar field Xi efficient preheating becomes possible even for rather small |£|. 

This paper is organized as follows. In the next section, we briefly review our model of R n inflation with the non- 
minimally coupled scalar field x- I n Sec. Ill, the preheating in the R 2 inflation model is investigated both by analytic 
approach and numerical integration. We show that analytical estimation based on the Mathieu equation coincides well 
with numerical results. In Sec. IV, the i? 4 inflation case is briefly discussed. We give our conclusions and discussions 
in the final section. 

II. BASIC EQUATIONS 

We propose a model with a higher-curvature term and a non-minimally coupled scalar field %, 

(2.1) 



-^R + a n R n - ^R X 2 - \(V X ) 2 - \™\ 



where k 2 /8ir = G = M pL is Newton's gravitational constant, a n and £ are coupling constants, and m x is a mass of the 
X field. Although this system contains only one scalar field x, it is reduced to a system described by Einstein-Hilbert 



action with two scalar fields by making conformal transformation to the Einstein frame 30 1 . We make the conformal 
transformation as follows: 

9i*v = e 2w sv, (2.2) 



where 



w = -4=^ = Jin [l + 2K 2 na n R n - 1 - ^V] 
V6 2 



(2.3) 



Then we obtain the following equivalent Lagrangian: 
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where the potential U n (4>, \) is defined by 



U n (cf > , X ) = e- 1 ^^ 
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±^R-a n R n + ^R X 2 + ±ml X i 



(2.5) 



Eliminating the Ricci scalar term by using Eq. (2.2), U n (<f>,x) can be rewritten as 
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where 



A„ = (n- 1) 



(2.7) 
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The scalar field ^ plays the roll of an inflaton field and the coupling constant a n (or A„) is constrained by (i) e- folding 
by the inflation and (ii) the density perturbation produced in the inflationary stage as follows. 

Since we need not take the x held into account in an inflationary stage, the potential U n can be approximated as, 



V n {cj>) ^ f«(<M) = A„e"^^ (e^ - l) "- 1 



(2.8) 



The form of the potential is quite different between the n = 2 case and others. Since the potential has a plateau 
for cj) > Mpl in the n — 2 case (See Fig. f), the inflaton rolls down very slowly on this plateau and can give the 
efficient inflation without making fine tuning on the initial condition of the inflaton field. This is one of the important 
properties of the R 2 inflation model. In the case of n > 2, however, V n ((f>) has a local maximum at 
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(2.9) 



and inflation becomes difficult to realize without fine tuning of the initial value because the inflaton easily falls down 
to the minimum of the potential. Let us investigate the initial and final value of <fi in the inflationary stage. First, 
the slow-roll parameters are described by the potential as 
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The inflationary stage ends when maxje!, |e 2 |} becomes of order unity. One can easily find that t\ increases to unity 
faster than |ea | , and the epoch when inflation ends can be estimated as 
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For example, 4> s = 0.188A/ PL w M PL /5 for n = 2, and 0/ = 0.108Af PL « M PL /9 for n = 4. Next, let us consider 
the value of <fi (= <fii) at the epoch of horizon exit when physical scales crossed outside the Hubble radius 60 e-folding 
before the end of inflation. By calculating e-folding number 



N= -k z 
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(2.13) 



and setting N — 60, we obtain the value of </ij as 
1.08A/ PL , 
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Mpl- (n > 4) 



(2.14) 



Note that <fii must be extremely close to the local maximum value 0* for n > 2. 

Next, we calculate the density perturbation produced in the inflationary stage and give the constraint of a n by 
comparing it with COBE data. The amplitude of density perturbation can be calculated as [[U[ 
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The COBE data requires 8h < 10 5 , hence a n is constrained to be 



a 2 > 10 9 , (n = 2) 

"„ = ln5 H 2{n - 1] > lO 1 ^"- 1 )^, (n > 4) 



(2.16) 



where 
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Note that a n needs to take a very large value especially in the n > 4 case. This is because the inflation field cj> must 
be very near to the local maximum of its potential initially in order to cause the sufficient inflation. For example, in 



the case of n = 4, the dimcnsionless coupling constant fa = a^M^ is restricted to be 



04 > 10 1 



(2.18) 



As the increase of n, we have to choose larger values of the coupling constant to fit the COBE data. However, since 
we are interested in preheating after inflation, it is worth investigating how x-particles are amplified apart from the 
constraint of the coupling constant even for the n > 4 case. 

After the inflationary stage ended, the universe enters the oscillating phase of the cj> field and x _ P ar ticles can be 
created because the \ field is non-minimally coupled with the spacetime curvature which is oscillated by the <f> field. 
Let us obtain the basic equations of the R n inflation model in the preheating stage. Since the x field is treated as a 
quantum field on the classical background of spacetime and the <j> field, the conformal factor in Eq. ( |2.2| ) includes a 
quantum variable. However, in order not to discuss quantum gravity, this term should be replaced with an expectation 
value (x 2 ), which may correspond to the number density of the ^-particle. Then we regard the conformal factor fl 2 
as 1 — r), where ry = £k 2 (x 2 ) fl32[ . 

Since we assume that the spacetime and the inflaton field are spatially homogeneous, we adopt the flat Friedmann- 
Robertson- Walker metric as the background spacetime; 



ds 2 = -dt 2 + a 2 {t)dx. 2 



(2.19) 



Hereafter we drop a caret since we argue only in the Einstein frame. From the Lagrangian (2.4), the evolution of the 
scale factor is described by 
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where a dot denotes a derivative with respect to time coordinate t. 
The evolution of the homogeneous inflaton field yields 



1+r? 



2 m x' 



;<x 2 



(2.20) 



3a 

a 



2y/Q 



>-K<f> 



1 Ve 



k4> 



((Vx) 2 



(e^^ -1 + 7?) 



+ A„e^-1 



+ ( n-2 



2(n-l) e ^- 1 + n + i W ^ 2) 



1 



= 0. 



(2.21) 



Since we postulate that the (j> field is spatially homogeneous, the fluctuation 5(f) is not considered. In the present 
model, the <f> field is a product of purely gravitational origin. Hence considering the fluctuation 8<j> is equivalent to 
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taking into account the perturbation of the metric. Recently, Bassett et al. |39 investigated the evolution of metric 
perturbation in a two-field model of a massive inflaton and a massless x. held interacting with the inflaton, and found 
that metric fluctuation is resonantly amplified in preheating phase. There will be a possibility that 5<j) grows by a 
parametric resonance also in the present model. Moreover, the rescattering effect of the x-particle and the 8<fi field 
are expected to make the x field be amplified. Although we do not consider the fluctuation in this paper, we should 
take into account this effect for a complete study of preheating. It is under consideration. 
The equation of the x field is expressed as 
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(2.22) 



where an index with a roman character denotes space coordinates. In order to study a quantum particle creation of 
the x field, we make the following mean field approximation with respect to x, which provides us a linearized equation 
for the quantum field x : 
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where V x = 2y / (x 2 ) d/d(x 2 )- Expanding the scalar field x as 
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(2.24) 



where a^ and al are the annihilation and creation operators, respectively, Xk obeys the following equation of motion: 
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Xk = 0. 



(2.25) 



The expectation value of x 2 is expressed with Xk as 



(X 2 ) = ^ J k 2 \ X k\ 2 dk. 



(2.26) 



Let us obtain the approximate equations of the <p an d X fields in order to estimate the preheating phenomenon 
analytically. First we rewrite Eqs. (2.2C) and ( [2.21 ) by following two approximations: 



I. The value of <fi is much smaller than Mpl, namely 



<M PL . 



(2.27) 



2. The backreaction effect by the x-particles is negligible, namely 



77 «f and ((V X ) 2 )«M| L . 
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Then, the evolution equation of the scale factor ( |2.20| ) is written by 
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where 
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(2.30) 



is a constant. On the other hand, the cj> field equation (2.21) can be written by 
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K n (j)— = 0. 



(2.31) 



Note that the potential of the <j> field in reheating phase is represented as V n {4>) « K n <fi~n=i . For n = 2, the potential 
becomes V2(4>) ~ K2<f> 2 , and the (f> field behaves as a massive inflaton. For n > 4, the potential is not regarded as 
the mass term and the curvature of the potential around its bottom becomes larger than that in the n = 2 case (See 
Fig 1). 

Making use of a time averaged relation: ((f) 2 ) = zrtr{V n (<f))) , we can easily find the time evolution of the scale factor 
from the relations (2.29) and (|2.3l|) as 



where to is the initial time of the coherent oscillation. Then, introducing a new time variable 
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and a new scalar field 
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the <fi field equation (2.31) is rewritten as 
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This equation is easy to be integrated and we find a conserved quantity, 
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Integrating Eq. ( 2.36 ) by the variable C = ("if) " l P 7rrT i T is represented as 
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For example, when n = 2, r is written as t — (l/v / 2~K^) cos" 1 C + ^o- Hence the oscillation of the ip field is sinusoidal: 
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For n > 2, however, the (p field does not oscillate sinusoidally but its behavior is rather complicated since the integral 



in Eq ( 2.37 ) is expressed by an elliptic or hyperelliptic integral. We investigate this issue for the n — 4 case in Sec. 
IV. 



Next we rewrite the x field equation (2.25) with the new time variable r. It yields 
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where 



Transforming the scalar field x^ as 
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Eq. ( |2.39D becomes 
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where the time dependent frequency u) 2 , is 
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y/6 dr 2 \ 3n~2 adr t/6 dr 
Since the frequency of the Xk field includes the inflaton field which oscillates coherently under two approximations 
( [2.27 ) and (2.28), the parametric resonance may be expected due to these terms even in the n > 4 case. We 



investigate x-particle production by the oscillating background field 4> especially when n = 2 and n — 4. We follow 
the semiclassical picture of preheating which is first considered in Ref. ]18[ , where the x~P ar ticles are created by 
quantum fluctuation initially, and they are treated as a classical field later as the x~P ar ticles are produced. We 
give an initial distribution of X^-state in our previous paper pq j. In the next section, we analyze the structure of 
preheating for the n = 2 case and compare it with numerical results. 

III. PREHEATING BY R 2 INFLATION MODEL 

In the n = 2 case, the reheating phase starts at 4> ss Mpl/5 and the (f> field begins to oscillate coherently around 
the minimum of its potential after <j) drops down to <\> < A:/pl/20. At this stage, as is found from Eqs. (2.32)-( 2.34| ), 



the scale factor evolves as a ss (t/to) 2 / 3 and r, tp are defined by r = t and ip — a 3 ^ 2 <f>. Since the if field oscillates 



sinusoidally as Eq. (2.38) with the mass 



-V^=-gp=, (3.1) 



the evolution of the <\> field for <\> < Mpl/20 is sinusoidal with decreasing amplitude 4> 



= $ sin mi, with $ = A Il L = M ^l_ , (3.2) 



which is equivalent to the case of the massive inflaton. Here we introduced the dimensionless time variable t = mi/27r, 



which represents the number of the oscillation of the inflaton field. Note that Eq. (3.2) is valid for t > 1. As the 
%-particles are produced, the backreaction effect by r/ and ((Vx) 2 ) is expected to appear. However, numerical 
calculations shows that rj and ((Vx) 2 )/»ri 2 AfpL do not exceed 0.1 even at the final stage of preheating, and the <p field 
oscillates almost coherently during the whole stage of preheating. 

Let us investigate the evolution of (x 2 ) both analytically and numerically. First, (x 2 ) is normalized as 

s2 > s w, = te) 2 i/ S2fo ' 2di = «^i/ p '«' 2 *' <"» 

where k = k/m and Xk — \frnxk- Hence (x 2 ) depends on the coupling constant a2- Here we adopt the lowest value 
«2 = 2x 10 9 which is determined by density perturbation. Note that in this case, m ~ 10 _6 Mpl ~ 10 13 GeV, which 



is almost the same value as that in the chaotic inflation model. Next, let us rewrite Eq. (2.42) to the form of Mathieu 
equation. Since the scale factor decreases so slow as compared with the time scale of the oscillation of the inflaton 
field for (/> < Mpl/20, we can approximate the frequency as 
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Here we did not neglect 77. Since 77 may increase up to 0.1, we should consider the effect from the backreaction terms 



£77 and 3ry in the X k field equation. Setting mt — 2z — 7r/2, Eq. (2.42) is reduced to the Mathieu equation: 
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77i x is normalized as 77i x = m x /m. Since we can find from Eq. (BJj) that A k always takes positive value even in 
the negative £ case [Remember 77 = £k 2 (x 2 )], the resonance bands are narrower than those in the model discussed 
in Ref. p6j when £ < case. However, since A k is not limited by A k > 2q, the resonance bands can be broader 
than those in the model of ^m 2 cf> 2 + \g 2 (j) 2 x 2 ■ Moreover q decreases as 1/t by the expansion of the universe. This is 
slower than in the massive inflaton plus coupled scalar field case, where q ~ 1/t 2 |26|]. This is due to the fact that the 



dominant oscillating term in Eq. (3.4) is not the </> 2 term but the cf> term, hence the period of the preheating becomes 
rather longer and we can expect an efficient resonance even if |£| is not very large. Below, we consider the case of 
£ > and £ < separately. 

A. Case of £ > 

1. massless x-particle case 

First, we investigate the production of the massless x-particlc for the case £ > 0. In the initial stage of preheating, 
we can neglect the £77 term. Hence A k is almost constant (A k ss 0) if we consider modes with small momentum k, 
which may be the dominant ingredients of the particle production independently on £. When the <fi field begins to 
oscillate coherently ((« 1), i.e., the analysis by using the Mathieu equation becomes valid, the initial value of q (= qi) 
is determined by £. After this, q decreases as 1/t by the expansion of the universe. Hence the Mathieu parameters 
trace the line near the g-axis horizontally on the Mathieu chart. When the chi field passes across the first instability 
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band, the resonance terminates. What determines the efficiency of the preheating is the initial value of q. When <& is 
large, the Floquet index is large and the growth rate of (x 2 ) is expected to be large. Moreover the period when the 
Mathieu parameters stay in the instability band becomes long, which means that the total abundance of (x 2 ) may be 
large. As a result, we can conclude that the large £ gives the efficient preheating as long as the backreaction effect 
can be neglected. Let us examine the evolution with the concrete values of £ numerically. 

For small £ (£ < 1), qi < 1. Hence, parametric resonance turns on in the narrow resonance regime and ^-particles 
are hardly produced because q decreases and the resonant range becomes narrower and narrower by the expansion of 
the universe. On the other hand, for £ > 1, the resonance can be expected. Let us consider several cases to compare 
the differences of the resonance structures. For £ = 3, qi ~ 3.61, and the resonance band is rather broad at the initial 
stage. In Fig. 3(a), we depict the evolution of (x 2 ) = (x 2 )/-^pl as a function of i for £ = 3. After i sw 1, (x 2 ) increases 
almost exponentially by the parametric resonance, and it reaches its maximum value (x 2 )/ = 2.21 x 10~ 8 when 
if = 5.61. [We summarize the data of numerical calculations in Table I.] After that, (X 2 ) decreases monotonically 
by the adiabatic damping of the Hubble expansion. Since the total amount of created particles is not so large, £77 is 
much smaller than unity during the whole stage and we can neglect the backreaction effect. In this case, the factor 
which terminates the resonance is the expansion of the universe. Ak and q both decrease and the resonance stops 
when q drops down to qf = 0.64 (if — 5.61). This evolution is illustrated by the solid curve in Fig. 2. For £ = 5, 
qi ~ 6.15, and more efficient resonance can occur. Actually, the final value of (X 2 )/ = 1-74 x 10~ 4 (when if = 6.27 



and qf = 0.98) is much larger than the case of £ = 3 (Fig. 3(b)). Since the term 36^rj in Eq. (3.6) is about 3.94 
at if = 6.27. Hence we can not neglect the backreaction in this case. The backreaction mainly affects the behavior 
of Ak- It makes Ak increase rapidly independently on the momentum k. Hence the Mathieu parameter deviates 
from the first instability band and the resonant band becomes very narrow. By this effect the parametric resonance 
terminates. For £ = 5, however, this deviation occurs when q drops down to about unity, hence we can say that the 
final value of (X 2 ) is mainly determined by the decrease of q. The case of £ rs 5 is the marginal case, whether the 
value of (x 2 )/ is determined by the decrease of q due to the expansion of the universe, or by the increase of Ak due 
to the backreaction. 

We can find that (x 2 ) increases with the same period as the oscillation of the inflaton field. This is different from 
the resonance in the chaotic inflation model with the massive inflation, where the (x 2 ) increases twice during one 
oscillation of the inflation p2| . In the chaotic inflation model, the Mathieu parameters are constrained as Ak > 2q, 
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and the square of the time dependent frequency of the x held is always positive. And the resonance occurs when the 
non-adiabatic condition: ui 2 <C dujk/ dt is satisfied. Since this condition is realized when <f> m 0, the resonance occurs 
twice during each oscillation of the inflaton field. In the present case, however, ui\ is not constrained to be positive 
but can become a negative value. In this case x-particles are produced mainly by negative coupling instability p5[ . 
As we can see from Eq ( |3.4D , u>\ becomes negative once for each oscillation of <f>. As a result, (x 2 ) increases once 
during each oscillation of the (f> field. This picture is independent on the sign of the coupling constant £. 

For £ = 10, q t ~ 12.52, and the final values are t f = 4.28, q f = 2.71, and (x 2 )/ = 9.95 x 1(T 5 (Fig. 3(c)). In this 
case, the value of (x 2 )/ decreases slightly compared with the £ = 5 case. As the x-particles are created, Ak increases 
because of the £77 term and the x field deviates from instability bands before q drops down under unity. This behavior 
is illustrated by the dashed curve in Fig. 2. Namely, in this case, the created x-particles themselves terminate the 
resonance. For £ > 5, the mechanism to stop the resonance is the same as in the £ = 10 case. One can estimate 
the final abundance (x 2 )/ by making use of this property. The maximal momentum to cause the negative coupling 



instability is given by Eq (3.4) as 



k 2 2 

max 

a 2 3ni 



3£(1 — 377) — m 



2 



9£7;-m 2 (3.8) 



x 2 

where k max = k max /m and fh x = m x /m. The relation ( p.q ) suggests that k'^ lax /a 2 decreases and approaches zero 
as the preheating proceeds. Actually, numerical calculation shows that K max ja 2 reaches zero at if in the case of 



£ > 5. For this reason, it is enough to set k 2 /a 2 = in Eq. (3J3) in order to obtain the analytical value of (x 2 )/- F° r 



£ > 5, parametric resonance terminates when Ak begins to increase rapidly as the increase of £77. Since the resonance 
bands are very narrow for the region Ak > 3<7, we adopt the criterion: When the Mathieu parameters pass across the 
Ak = 3q line, the preheating terminates. Then, we find the final abundance of the x-p&rticles as 

Note that (x 2 )/ includes two variables: £ and if. Generally, as £ increases (£ > 5), if decreases because the production 
rate of x-particles becomes large and the x field deviates from the instability bands in the early stage. Numerically, 
(x 2 )/ is the slowly decreasing function of £ for £ > 5 and this suggests that the £ effect plays a more important 



role. The analytic estimation by Eq. (3.9) shows good agreement with numerical calculation. For example, when 
£ = 10, the estimated value is (x 2 )/ = 9.02 x 10~ 5 and the numerical value is (x 2 )/ = 9.95 x 10~ 5 . For £ = 100, 
the analytic value is (x 2 )/ = 2.87 x 10~ 5 , since if = 1.15. This is close to the numerical value (x 2 )/ = 2.51 x 10 -5 . 
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This coincidence shows the validity of the criterion and the picture which we adopted. For £ > 100, If approaches a 
constant value if ~ 1, and (x 2 ) f decreases as l/£. Other examples are listed in Table I. We also show the final value 
of (x 2 )/ as a function of £ in Fig. 4. For 1 < £ < 5, (x 2 )/ is the increasing function of £ and is determined by the 
decrease of q due to the Hubble expansion. On the other hand, (x 2 )f is a slowly decreasing function of £ for £ > 5, 
where the backreaction effect becomes important. 

We can conclude that (x 2 ) f takes maximal value (x 2 )max — 1-74 x 10~ 4 when £ « 5 for £ > for the massless 
X-particle case. 

2. massive x-par£ic/e case 

Whether the x-P ar ticl e with a mass heavier than the GUT scale can be produced or not is one of the most 
important issues of the preheating, which is related to GUT scale baryogenesis. One can easily obtain the properties 



of the preheating by the massive x field and difference from the massless case from Eqs. (3.6) and Eq. (3.7). Since 
the x - field couples with the inflaton field by £ and m x terms as we can see from Eq. ( |2.25| ), one may consider that 
parametric resonance would occur by only the mass term without the non-minimal coupling. This, however, is not 
the case. In the case of £ = 0, Ak > 4m 2 and |cj»j| « 0.4 (m 2 + i). This means that Ak is at least 10 times larger than 
q when i ss 1, in which case the resonance bands are very narrow. Moreover, even if the x field stays in an instability 
band initially, Ak and q decrease by the expansion of the universe and soon deviate from the instability bands. As a 
result, it is difficult to produce the massive x _ P ar ticle enough without the non-minimal coupling. We confirmed this 
by numerical calculation. 

For the m x ^ case, the existence of the m 2 term constrains Ak > 4m 2 and decreases the value of q, hence the 
resonance is suppressed as compared with the massless case. First, consider the case of m x — m (m x — I). Ak takes 
the value more than 4 due to the mass effect, and the modes with the momentum which contribute to the resonance 
mainly exist in the second instability band of the Mathieu chart (See Fig. 2). Since there is no contribution from the 
first resonance band, which has larger Floquet index than the second one with the same q value, the larger value of £ 
is expected to be needed for the ^-particle production to occur. Numerically, when £ < 3, the parametric resonance 
does not come about. In the case of £ = 5, the resonance occurs, but this process is very weak (Fig. 5(a)). For 
t > 1, (x 2 ) increases, but (x 2 )f is smaller than the initial value. When £ = 10, the parametric resonance after i s=s 1 
is rather effective (Fig. 5(b)). The final values are if = 9.22, q f = 1.31, and (x 2 )/ = 1.17 x 10~ 5 (See Table II). 
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Since the contribution from the backreaction term 36£?y in Eq. (3.6) takes 1.06 as a final value, the deviation from the 
instability band by the increase of Ak begins to be important. When £ = 20, the final values are if = 3.23, q/ = 7.07, 
and (x 2 )/ = 5.19 x 10~ 5 (Fig. 5(c)). In this case, the backreaction effect by the £77 term completely determines the 
final x-particle abundance because qf = 7.07 when resonance ends. Numerically, (x 2 )/ takes the maximal value of 
(x 2 }max ~ 5 x 10~ 5 for £ = 20 ~ 30. For £ > 30, (x 2 )/ is a slowly decreasing function of £. Generally, when the £77 
suppression effect plays the relevant role, one can estimate (x 2 )/ in the same way as the massless x _ P ar ticle case. By 
adopting the criterion: Ak = 3g when the resonance terminates and setting k 2 /a 2 = in Eq. (3J3), we obtain 



(x 2 )/ 



72< 2 



6£-l 



K 



(3.10) 



.2(7rt> + 1) 

Apparently, the mass effect suppresses the final x _ P ar ticle abundance. Let us compare this estimation with numerical 
results. For example, when fh x = 1 and £ = 20, the numerical value is (x 2 )/ = 5.19 x 10~ 5 (at if = 3.23) and the 
analytic one is (x 2 )/ = 4.80 x 10~ 5 ; for £ = 100, the numerical one is (x 2 )/ = 2.78 x 10~ 5 (at if = 1.18) and the 
analytic one is (x 2 )/ = 2.76 x 10~ 5 . When fh x — 1, the mass term in Eq. ( |3.10] ) can not be neglected as compared 
with the former term for the £ < 50 case. In this case, the mass effect is one of the important factors to determine 
the final abundance of the \ particles. On the other hand, for £ > 50, the mass effect can be neglected and does not 
play a relevant role for the final abundance, although the growth rate and other properties during the resonance are 
different from the massless case. In Table II, we show the analytically estimated and numerical values of (x 2 )/ in 
various parameters. We can see that the analytic results give good agreement with the numerical ones. 

For the large £ where the backreaction becomes important, the right hand side (r.h.s.) in Eq. ( |3 .10 ) must be a 
positive value, hence fh x is constrained to be 



* - 2(ntf + 1) 



K % o/r , V (3-H) 



Since tf approaches to the constant value of tf fa 1 as £ increases, we can obtain the constraint of the mass of the 
created x _ particle for large £ as 



'»> 



< v 7 ?- (3.12) 



From this constraint we find that if £ is more than 100, x - P ar ticle whose mass is of order m x = 10m ~ 10 GeV 
can be produced by the non-minimal coupling in the R 2 inflation model. In order to produce the GUT scale gauge 
boson whose mass is about m x = 10 3 m ~ 10 16 GeV, £ is needed as at least 10 6 . Namely, unless £ takes a very large 
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value, the GUT scale boson is hard to be created. Moreover, even if £ takes a very large value, since the Mathieu 
parameters Ak and q deviate from instability band due to mass effect from the very beginning, the final abundance of 
the x-particles must be very small as estimated by Eq. (3.10|) as (x 2 )/ « !/?■ Finally, we show the value of (x 2 )/ as 



a function of m x in the case of £ = 100 in Fig. 6. (x 2 )/ is the quadratically decreasing function of m x as is expected 
from Eq. (3.10p. We find that the x-particle whose mass is more than m x > 10m can not be produced, justifying the 



estimation (3.12|) 



B. Case of £ < 

1. massless x-particle case 

Let us investigate the £ < case for the massless x-particle. In the £ < case, q takes negative value as is found 



by Eq. (3/7). Since the Mathieu chart is symmetrical with respect to A^-axis, it is enough to consider the absolute 
value of q: 

4 



191 = 3^" 



3|£l(l-3r/) + i 



(3.13) 



Note that \q\ is larger than in the positive £ case because of the last term +1/2, which makes the resonant more 
efficient. Even if £ = —1 (g, fw 1-5), parametric resonance evidently occurs, although the final abundance is very 
small: (X 2 )/ = 2.63 x lO" 11 (see Fig. 7(a)). For £ = -3 (% « 4.0), the final values are if = 7.02, q f = 0.57 and 
(X 2 )/ = 8.71 x 10~ 5 . Hence parametric resonance is much more efficient than in the £, — 3 case. One characteristic 
property in the £ < case is that (X 2 ) can increase before the <p field enters the coherent oscillation stage. This is 
because the time dependent frequency u>\ of Eq. (p. 43) takes negative value initially for k — 0, and (x 2 ) grows due 



to negative coupling instability. For example, in the £ = —3 case, (x 2 ) decreases slightly initially by the expansion of 
the universe, but begins to increase after that (Fig. 7(b)). This is another factor where the negative £ case gives more 



efficient preheating than the positive one. After {wl, the approximation (2.27) is valid and the analytical estimation 
by the Mathieu equation is relevant. The final value qf is small for £ > —3, we can neglect the back reaction effect 
and the resonance terminates by passing across the first resonant band. For £ = — 5, however, the final values are 
if = 3.80, qf — 1.57 and (x 2 )/ = 2.63 x 10~ 4 (Fig. 7(c)). Since qf is larger than unity, the backreaction by the £?y 
term in Eq. (|3.6|) becomes important and one can estimate the final abundance (x 2 )/ for £ < — 4 by the same method 



as the £ > case. By using the criterion Ak = 3\q\ for Eqs. (3.6) and (3.13), (x )/ is estimated as 
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(x 2 )/ 



247^1(7^-1) 



m 



(3.14) 



This shows that the maximal value of (x 2 )f becomes larger than that of the £ > case in Eq. (3/7). Actually, the 
value of (x 2 )f in the range —15 < £ < — 4 is larger than the maximal value for £ > case (See Fig. 4). For £ = —10, 



(X 2 )/ estimated by Eq. ( |3.14j ) is (x 2 )/ = 1.74 x 10~ 4 (at t f = 2.78) and the numerical value is (x 2 )/ = 1.78 x 1(T 



The analytic estimation (3.14) gives quite good agreement with the numerical results for —20 < £ < — 4 (Sec Table 
III). 



For £ < — 20, however, the final point tf becomes less than unity, and the estimation (3.14) based on the Mathieu 



equation is broken. In this case, since tof, takes negative value for most k modes, (x 2 ) and \r]\ increase rapidly. Then, 



the term e 3 K< ^ — 1 4- 377 in Eq. (2.43|) decreases and becomes negative while lo\ becomes positive for any <fi, which 



means the negative coupling instability does not occur any more. Note that the field equation ( 2.21 ) includes a 



;T, 



similar term : e~3~ re * — 1 + 77. Although the r\ term is expected to work as the backreaction on the 4> field, the effect 
on the x field equation is dominant due to the difference by factor 3. After the inflaton field enters the coherent 
oscillation stage, A^ is much larger than \q\ and the further resonance does not come about. By making use of this 
picture, we can estimate the final abundance of (x 2 )/ f° r £ < — 20 as 

where <\>f = </>(£/). For example, when £ = —50, numerical values are if = 0.58, (x 2 )/ = 2.82 x 10 -4 and <f>f — 



0.128AfpL- Then, (x 2 )/ estimated by Eq. (3.15) is (x 2 )/ = 1-83 x 10~ 4 . In spite of the naive estimation, it gives a 
value close to the numerical one. We show the final abundance (x 2 )/ for £ < in Fig. 4. We can find that (x 2 )/ takes 
its maximum value (x 2 )max ~ 3 x 10~ 4 for £ f=a —4 (and £ w —40), which is larger than £ > case. For £ < — 40, 
(x 2 )/ begins to decrease by the factor: 1/|£|. 

2. massive x-po-rticle case 



We investigate the £ < 0, m x ^ case as the final case in this section. Since A^ and \q\ are expressed as Eq. ( |3.6| ) 
and 



37rt 



3|£|(l-3r/) + m : 



., 1 

+ - 

* 2 



(3.16) 



the mass effect makes both Ak and \q\ larger than the massless case. Although the large q value gives an efficient 
preheating, the increase of Ak is more significant than the increase of |g|, and x _ particle production is suppressed by 
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taking the mass effect into account. For example, when rh x = 1, the resonance of x-P ar ticles can not be seen for 
— 3 < £ < 0. Consider the case of £ = — 5 and fh x — 1 (See Fig. 8(a)). Although the growth rate is small, (x 2 ) is 
slowly increasing up to (x 2 )f = 3.24 x 1CP 8 (at if — 7.74). In this case, since g/ = 0.90 and 36£?7 <C 1, the resonance 
terminates by the expansion of the universe. When £ = — 10 and rh x = 1, the resonance is efficient and the final 
values are (x 2 )f — 1.13 x 10~ 4 , if — 2.69 and qf — 4.61 (Fig. 8(b)). In this case, the increase of Ak due to x-particle 
production determines the final abundance. By using the same criterion as in the previous cases, we can estimate the 
value of (x 2 )f as 



(x 2 : 



i 



/ 



72tt£ 2 



6|£|+1 



2(ni f - 1) 



(3.17) 



In the £ = —10 case, the estimated value is (x 2 )f = 1-37 x 10~ 4 , which is close to the numerical value. Eq. (3.17) is 
relevant for —30 < £ < — 10 because if is larger than unity. However, for £ < — 30, where the resonance terminates 
before the inflaton field begins to oscillate coherently, we should use the estimation 

1 



(x 2 ; 



■m \ 



f 



e V3**/ - 1 * 

24^i\ ml' 



(3.18) 



which can be obtained similarly to Eq. ( 3.15| ). Note that Eq. ( 3.18|) is valid when m 2 is not as large as |f |. If m 2 is 



large as compared with |£|, uf. becomes positive at the initial stage by Eq. ( 2.43 ) and can not give the negative coupling 
instability. For example, when £ = —50 and fh x — 1, the numerical values are (x 2 )f = 2.79 x 10~ 4 , <f>f — 0.128Mpl 



at tf — 0.58, while the estimated value by Eq. ( 3.18] ) is (x 2 ) f — 1-81 x 10 (See Table IV for other examples). In 



this case, as is found by Eq. ( 3.18 ), the mass term does not affect a significant role. However, if m x is large, (x 2 }f 
is significantly suppressed by a mass effect (Compare with £ = —50, rh x — 7 case). Since x-particle production is 



possible naively when the r.h.s. in Eq. (3.18) takes positive value, the produced mass of the x-p&rticle is bounded as 



m x £ 



^(eV?"*/-!). 



(3.19) 



4>f is smaller than the initial value Mpl/5 for large £. This means that m x is constrained to be 



™* < 2v/j|f. 



(3.20) 



The massive x _ particle whose mass is of the order m x ~ 10 GeV can be generated even if £ ~ —5, but with the 
increase of m x , the x-particle is hard to produce because m x is strongly restricted by Eq. ( 3. 20| ). In particular, |£| 
still needs to take a very large value |£| > 10 5 in order to produce the GUT scale gauge boson: m x ~ 10 16 GeV. As a 
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result, it seems difficult to produce the GUT scale particle by parametric resonance in the R 2 inflation model in both 
the £ > and £ < cases. 

IV. PREHEATING BY R 4 INFLATION MODEL 



In this section, we consider the R A inflation model as the next example. In the R 4 inflation model, inflation is 
hardly realized because a fine tuning of the initial value of the inflaton field is needed and the coupling constant 0:4 
is constrained to be very large |34j. It is worth, however, investigating whether a preheating stage exists or not and 
how the preheating process proceeds if it occurs. In this case, the value of the inflaton when the reheating process 
turns on is « Mpl/9. After <f> decreases under <f> « Mpl/20, the scale factor evolves asa« (t/to) , where r and 
ip are defined as r = 2£q/3 (i/io) an d </? = a 9 ^ 5 4>, respectively by Eqs. ( 2.32 )-( 2.34 ). As for the ip field, integrating 



the r.h.s. in Eq. ( 2.37 ), we find that r is expressed by the combination of the complete elliptic integrals of the first 
kind F (9, 1/V2) and the second kind E (9, I/a/2) §^: 



ne ^)- 2E {"-T2 



±- 



2 (Ki_ 

E 



3\V4 



( T - To), 



(4.1) 



where 



1/4 



(4.2) 



Since F (9, 1/V2) and E (9, 1/^2) are well approximated as F (9, l/%/2) « 1.18035 and E (9, 1/^2) sw 0.85986*, we 
obtain the evolution of the ip field as 



9' : 



<Pq cos 



L Vo 



(4.3) 



where c = 1.2360 is a constant and po = (E/K4) 4 represents the amplitude of the p field. Note that the variable 



TO4 



1/3 



(4.4) 



plays the role of the "inflaton mass" for the R 4 inflation model when the inflaton field oscillates coherently. 
Next, consider the evolution of the x field. The time dependent frequency in the present model is expressed as 



-A 



£n 2 (e^ 



- 1 + ?/ 



vfi 



^ - 1 + -ri + m 



k 2 f d(f> 



K d 2 



y/6dr 2 6 \dr 



(4.5) 



where we used the approximation n-£p ^> 3^. Further, adopting the approximations: Eqs. ( 2.27 ) and ( 2.2§| ), the 
frequency can be written by 
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UJ l~ a § (^2 + m x) + (3) a H^- I) nK 4 (fil cos [m 4 (T-T )}, (4.6) 



where we used Eq. (4.3). Setting 2z = m^r — t ) + 7r, Eq. (2.42) is reduced to the Mathieu equation (3.5), with 

'P 



A k=^{y+<)' ( 4 - 7 ) 

q = c / U-^ja- C icp , (4.8) 

where c' = fi\/¥ = 21.4338, fc = k/m^, fh x = m x jra^ and (^0 = Vo/^pl- Both of the Mathieu parameters decrease 
by the expansion of the universe in the R model. Note that this relation can only be used for 77 <C 1. As x-particles 
are produced, we have to take into account the suppression effect by the increase of \rj\. As the R 2 inflation model, 



the difference between the £ > and £ < cases appears in the sign of the term 1/6 in Eq. ( |4.8| ). This difference, 
however, is not important for |£| 3> 1. 

The expectation value (x 2 ) is normalized as 

<X 2 >^=feVi fk 2 \x k \ 2 dk= — -4tT7T ^L [k 2 \x k \ 2 dk, (4-9) 



M| L UfpiJ 27r 2 7 IAK| 16(127T^ )^M 2tt 2 



where x& = ^Jm^Xk- Since /3 is strongly constrained as Eq. (2.18), the normalized value of (x 2 ) is very small initially. 
However, since (x 2 ) can grow exponentially by the parametric resonance, we can expect a considerable amount of the 
X-particles at the final stage of the preheating. We use a dimensionless time parameter defined by t = Mp^(3~et in 



numerical calculations. Counting time from <j) ss Mpl/9, the evolution of the ip field is described by Eq. (|4.3| ) with 



</>o w 0.05M PL (at t w 20). Then Eq. ( f4.8| ) is rewritten as 



.»u>™(< -£)(£)"*. (4.10) 

We find from this relation that q decreases faster as compared with the R 2 inflation model. 

First, consider the £ > case for the massless %-particles. For < £ < 3, (x 2 ) does not increase because the initial 
value of q is small. For £ > 5, however, we can see that (x 2 ) grows by the parametric resonance. Consider the case of 
£ = 10 (See Fig. 9(a)). Numerically, the final values are if = 161 and (x 2 )/ = 5.89 x 10~ 27 (See Table IV). Although 
the total amount of created x _ P ar ticles is not large because of the small initial value, parametric resonance evidently 
occurs due to the oscillating ip field. In this case, since rj is much smaller than unity, the backreaction effect due to 
X-particle production does not affect the x _ field evolution. When the Mathieu parameters with low momentum, which 
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contribute the resonance mainly, pass across the first instability band by Hubble expansion, the resonance terminates. 
In the case of £ = 30, the final values are if = 185 and (x 2 ) — 1 -29 x I0 -5 (Fig. 9(b)). In spite of the small initial 
fluctuation, (x 2 ) increases up to (x 2 )/ ~ 10~ 5 by copious production of the x-particles. In the case of £ = 50, the 
final values are if = 140 and (x 2 )/ = 8.28 x 10~ 6 (Fig. 9(c)). Although the growth rate of (x 2 ) becomes larger than 
the £ = 30 case because of the large initial value of q, the final abundance is almost the same. This means that the 
suppression effect by the backreaction is significant when the fluctuation increases up to (x 2 ) ~ 10~ 5 . In such cases, 



we can estimate the final abundance of x _ P ar ticles as follows. When i)<l, the time dependent frequency (4J3) can 



become negative for 4> < and the x - P&rticle can be produced by the negative coupling instability. However as the 



X-particles are produced enough, the second term in the square bracket of Eq. (4.5) becomes positive for any value of 



We can roughly estimate that the resonance finally terminates when 77 grows up to 

where we used the relation $y ss f -£- \ <p$. $ represents the amplitude of 0. Namely, 

2, . 1 /3/t/ 



^^lolVnt' *°- (4 ' 12) 



Let us adopt the values of ipo = 0.05 and to = 20. For the £ = 30 and £ = 50 cases, (x 2 )/ estimated by Eq. (4.12) 
are (x 2 )f = 5.79 x 10 -6 and (x 2 )/ = 5.28 x 10~ 6 , respectively (Other examples are listed in Table V). In spite of the 



naive estimation of Eq. (4.12|), it gives fairly close values to the numerical ones. We did not take into account the 



d 2 (j)/dT 2 term in Eq ( |4.5| ) for the estimation. This is one of the reasons why the analytical values are smaller than the 
numerical values. Note that the r.h.s in Eq. (1.12|) can be regarded as the function of two variables £ and tf. With 



the increase of £, numerical calculation shows that tf decreases. For 25 < £ < 100, the decreasing rate l/£ of (x 2 )f 

almost balances with the increasing rate (tf/to) 2 ■ Hence (x 2 ) f takes almost constant value (x 2 )/ ~ (1 ~ 3) x 10~ 5 

(See Fig. 10). To be more precise, (x 2 )/ takes the maximal value (x 2 )/ ~ 3 x 10~ 5 for £ ss 40. For £ > 100, the 

_ 3 
decreasing rate l/£ surpasses the increasing rate ( j- ) and (x 2 )/ decreases with the increase of £. 

Next, we consider the £ < case for the massless x-particles. In this case, since the absolute value of q is represented 



as 



W ="'! !£! + £)(£) 2 ^o, (4.13) 



the initial value of \q\ is slightly larger than the £ > case. Moreover, since u 2 in Eq. (4.5) can take negative value from 
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the first stage of preheating, (x 2 ) grows rather effectively. Numerically, while (x 2 ) does not increase for —3 < £ < 
at all, the growth of (% 2 ) is expected for £ < — 3. For example, in the case of £ = —5, (x 2 ) increases exponentially 
and reaches the final value (x 2 )/ = 1-95 x 10~ 33 at if = 100 (Fig. 11(a)). At the initial stage, (x 2 ) decreases slightly 
by the Hubble expansion, and then increases due to negative instability from t ss 10, which is earlier than the £ > 
case. The smaller the value of £ we adopt, the larger the final abundance (x 2 )/ becomes, and it reaches a plateau 
when £ « -20. In the case of £ = -20, the final values are if = 156 and ( X 2 )f = 3.08 x 10" 5 (Fig. 11(b)). As £ 
further decreases under £ w —20, the suppression effect by ^-particle production begins to work, and (x 2 )/ does not 
increase while the growth rate becomes large as the decrease of £. Numerically, (x 2 )/ takes almost constant value 
(x 2 )/ fa (2 ~ 5) x 10~ 5 for the range —100 < £ < — 20. For example, in the case of £ = —50, the final values are 
if = 82 and (x 2 )/ = 2.18 x 10~ 5 (Fig. 11(c)). We can estimate (x 2 )/ in the same way as the £ > case, which yields 



When £ = —50, (x 2 )/ estimated by Eq. ( 4.14 ) is (x 2 )/ = 1-18 x 10 5 , which is close to the numerical value (x 2 



>/ 



2.18 x 10 5 . In the case of £ < — 100, the factor 1/|£| becomes dominant and (x 2 )/ decreases with the increase of |£|. 



In Table VI, we show the numerical and analytical values by Eq. ( 4.14 ) in various cases. We also depict the numerical 
value of (x 2 )/ as a function of £ in Fig. 10. We find that (x 2 )/ takes almost constant value for —100 < £ < — 20. 
The achieved maximal fluctuation in the R A inflation model is {x 2 )max ~ 5 x 10~ 5 at £ fa —35, which is smaller than 
that in the R 2 inflation model. 

Finally, we take the mass of the x-particles into account. In the R 2 inflation model, Ak is constrained as Aj~ > 4m 2 
throughout the preheating. In the R A inflation model, although A^ increases by the mass term similarly, it can 
decrease due to the prefactor a~ 6 ' 5 and the system might give a resonance which is not inefficient compared with 
the massless case. However, we have to note that q value decreases by the Hubble expansion simultaneously. If q 
decreases faster than A^, the parametric resonance is suppressed as in the R 2 inflation model. Unfortunately, this is 
the case. As one example, we show the (x 2 ) evolution in the case of £ = —50, m x = 7714 in Fig. 12. We find that the 



growth of (x 2 ) is strongly suppressed by the mass effect. Although Ak and q in Eq. (4/7) and (|4.8| ) are not written 
by the forms which include the term 77, we can roughly estimate the maximal value of m x by using the criterion: 
Ak < 3q. Considering the k = mode, fh x is constrained as 
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The above estimation gives a much lower upper bound than that in the R 2 inflation case. From Eq. ( Ill , the 



X-particle whose mass is of the order m x ~ 7714 ~ 3 x 10 22 A/pl can be produced if |£| > 100. We conclude that an 
unnaturally large value of |£| is needed even if GeV scale bosons are to be produced. 

V. CONCLUSIONS AND DISCUSSIONS 

In this paper, we analyze preheating with non-minimally coupled scalar field \ m the higher-curvature inflation 
model. We have examined properties of resonance, especially for R 2 - and _R 4 -inflation models. In the R 2 model, 
i? 2 -term provides us an effective inflaton field (ft, and inflation is realized through the flat plateau for (ft > Mp^. In the 
reheating phase, the (ft field behaves as a massive scalar field. Although the evolutions of scale factor and (ft field in the 
present model are almost the same as those for a massive inflaton with V{(ft) = ^m 2 (ft 2 in the non-minimal coupling 
■i£i?X 2 case, the structure of resonance is different. In the R 2 model, since the resonance parameter q decreases as 
1/t due to the expansion of the universe, which is slower than that in the model with V((ft, x) — 2 m2( ^ 2 + ^C-^X 2 : we 
do not need so large a value of |£| as |£| > 10 for an effective resonance. In the £ > case with a massless X"P ar ticle, 
if £ is larger than £ > 1, X"P ar ticle production is possible against the diluting effect by cosmic expansion. When 
1 < £ < 5, the total amount of created x-particles (x 2 )/ grows with the increase of £, because the initial value of q 
becomes larger. However, for £ > 5, (x 2 )/ slowly decreases as the increase of £ because of the suppression effect by 



X-particle production. Hence, in the positive £ case, (x 2 )/ takes the maximal value \/(x 2 ) ~ 1-5 x 10 17 GeV at 
£ ~ 5. When £ is negative, the resonance structure does not change so much except for the small increase of the \q\ 
value. In this case, the x-particle production becomes possible for £ < — 1. For —4 < £ < — 1, (x 2 )/ increases with 



the decrease of £, and (x 2 )/ takes maximal value \J {x 2 ) ~ 2 x 10 17 GeV at £ w —4. For the case of £ < — 20, a 
X-particle fluctuation rapidly grows with the passage of time and reaches the maximal value (x 2 )/ before the (ft- field 
begins to oscillate sinusoidally. This means that in the £ < case, the effective x-particle production occurs due to a 
negative coupling instability irrespective of the coherent oscillation of 0-field. Numerically, (% )/ slowly increases for 
—40 < £ < — 20 with the decrease of £, and the value of (x 2 )/ m the £ « —40 case is almost the same as the £ « —4 



case. In summary, in the R 2 model, the achieved maximal value of x-P ar ticle fluctuation is y/(x 2 } ~ 2 x 10 17 
GeV for £ w —4 (and £ « —40). As for the massive x-particle case, a parametric resonance is suppressed by the mass 
effect. We find the relation for the massive x-particle production to occur: m x < y/£ for £ > and fh x < 2y/\£\ for 
£ < 0. Hence it is difficult to produce GUT scale gauge bosons unless |£| is as enormously large as |£| > 10 5 . 
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As for the R A model, inflation is difficult to realize unless selecting the very large coupling constant a^ with 
spacetime curvature. Nevertheless, we have examined the resonance structure, because we are interested in whether 
a preheating stage generally exists or not for n > 2. Since the potential of the <\> field is written as V((f>) ~ 4>z for 
4> < Mpl/20, the oscillation of the (f> field is not exactly sinusoidal. However, we find that the behavior of cf> field 
is well approximated as a sinusoidal oscillation, and the equation of the \ field is reduced to the Mathieu equation. 
In the R 4 case, the resonance structure changes and q decreases as q ~ t~?, which is faster than the R 2 case. For 
£ > 5 and £ < — 3, the fluctuation of %-particles grows with the passage of time. However, (x 2 )/ is rather small 
for |£| < 20, because a parametric resonance is not so effective. For 20 < |£| < 100, (x 2 )/ takes the almost constant 
value \J{x 2 ) f ~ (4 ~ 8) x 10 16 GeV, and for |£| > 100, (x 2 )/ decreases with the increase of |£|. The maximum value 
achieved in the i? 4 model is \/(x 2 ) t ~ 8 x 10 16 GeV at £ sw —35. Taking into account the mass effect of x-particles, 
the final fluctuation is more strongly suppressed as compared with the R 2 inflation model. 



As for R n (n > 4) model, the coupling constant a n is constrained to be very large as Eq. (2.16). Although this 



model is rather unnatural, we should stress that it is still possible to investigate the resonance structure by the same 



method as in the n — 4 case. The oscillation of the <j) held around the potential Eq. (2.8) is approximately written 
by the sinusoidal form, so the fluctuation of ^-particles can grow by a parametric resonance. Hence, in this sense, a 
preheating stage generically exists in thw higher-curvature inflation model with a non-minimal coupling, irrespective 
of the problem of the large coupling constant. 

In a higher-curvature inflation model, inflation is realized through the coupling with spacetime curvature. Because 
of this, the inflaton field is the product of a purely gravitational source. Hence for the complete study of preheating, 
we should take into account the metric perturbation. It was pointed out in Ref. [|36| fllfl that the metric perturbation 
undergoes the effect of parametric resonance in a reheating phase. Although it is rather difficult to study preheating 
in the present model by taking into account metric perturbations, it is of interest how the preheating dynamics are 
modified by this effect. 

In our study of amplification of a scalar field x coupled non-minimally to a spacetime curvature R (^£i?x 2 ) m 
a higher-curvature inflation model (R + a n R n ), we have found that the value of £ does not need to be so large as 
in the case of massive inflaton plus non-minimally coupled scalar field x to have an effective resonance. In other 
gravitational theories such as the Brans-Dicke theory p3] and the induced gravity W3[ , we may have different types 
of preheating dynamics. Recently, Mazumdar and Mendes B3] studied preheating in the Brans-Dicke theory with a 
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model of massive infiaton plus a minimally coupled scalar field a and found that the Brans-Dicke field as well as the a 
field undergoes the effect of parametric resonance. It is of interest to investigate the amplification of a non-minimally 
coupled scalar field in Generalized Einstein Theories. Furthermore, although we mainly concentrate on the x-particle 
production in this paper, we may discuss an infiaton particle production with a model of some infiaton potential and 
non-minimally coupled scalar field <p as well. These issues are under consideration. 
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TABLE I. The final values (x 2 )f obtained by the analytical estimation and by the numerical calculation for £ > 0, m x = in 
the case of the R 2 inflation model. The analytical estimation gives a good approximation to the numerical results, (x 2 ) numerical 
takes the maximal value at £ ~ 5. We also show the time if when the resonance stops and the value of 36^7/ which indicates a 
suppression effect. For large values of f , the suppression effect by x~ particle production is significant. 



£ 


*/ 


\X /analytic 


\X /numerical 


^vZ^Tfn umerical 


3 


5.61 


— 


2.21 x 10" 8 


1.80 x 10~ 4 


5 


6.27 


1.24 x 10~ 4 


1.74 x 10~ 4 


3.94 


10 


4.28 


9.02 x 10~ 5 


9.95 x 10~ 5 


9.00 


30 


2.15 


5.67 x 10~ 5 


7.50 x 10~ 5 


6.10 x 10 1 


50 


1.30 


5.20 x 10~ 5 


2.90 x 10~ 5 


6.56 x 10 1 


100 


1.15 


2.87 x 10~ 5 


2.51 x 10~ 5 


2.27 x 10 2 


500 


1.06 


6.13 x 10~ 6 


7.96 x 10~ 6 


1.80 x 10 3 



TABLE II. The final values (x )f obtained by the analytical estimation and by the numerical calculation for £ > 0,m x 7^ 
in the case of the R 2 inflation model. The mass effect suppresses the final x _ P ar ticle abundance, and massive x _ P ar ti c ' e which 
does not satisfy the condition Eq. (3.14) can not be created. 
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\X /analytic 
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3.80 x 10 -12 
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1.17 x 10~ 5 
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4.80 x 10~ 5 


5.19 x 10~ 5 


3.23 
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2.80 x 10~ 5 


2.51 x 10~ 5 


2.08 


100 
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2.36 x 10~ 5 


2.32 x 10~ 5 


1.21 


500 


10 


3.87 x 10~ 6 


3.80 x 10~ 6 


1.18 



TABLE III. The final values {x 2 )f obtained by the analytical estimation and by the numerical calculation for £ < 0,m x = 
in the case of the R 2 inflation model. The analytical estimation gives a good approximation to the numerical results. We also 



show the time tf when the resonance stops and the value of 36^77 which indicates 
maximal value at £ w —4 and £ ~ —40. 



suppression effect, (x ) 
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1.83 x 10~ 4 


2.82 x 10~ 4 
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9.13 x 10~ 5 
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1.10 x 10 3 


-500 


0.37 


3.19 x 10~ 5 


5.31 x 10~ 5 


1.20 x 10 4 
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TABLE IV. The final values (x 2 )/ obtained by the analytical estimation and by the numerical calculation for £ < 0,m x 7^ 
in the case of the R 2 inflation model. The mass effect suppresses the final ^-particle abundance, and massive x _ P ar ticle which 
does not satisfy the condition Eq. (3.22) can not be created. 
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TABLE V. The final values (\ }/ obtained by the analytical estimation and by the numerical calculation for £ > 0,m x 
in the case of the R 4 inflation model. 
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TABLE VI. The final values (x 2 )/ obtained by the analytical estimation and by the numerical calculation for £ < 0,ra x 
in the case of the R 4 inflation model. 
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Figure Captions 



FIG. 1: The potential in the equivalent system to the R n gravity thoery. The solid and 

dotted curve denote the n — 2 and n — 4 cases respectively. The potential of 
the n = 2 case is flat enough for cj> > Mpl so that inflation occurs, but for the 
n > 2 case, the potential has a local maximum at Eq. (2.9) and inflation is 
rather difficult to realize. 

FIG. 2: The schematic diagram of the Mathieu chart and the typical paths for two types 

of resonance. The lined regions denote the instability bands. Generally, The 
growth of x fluctuation terminates in two distinct ways. In the first case, when 
|£| is not so large, resonance ceases by universe expansion before x _ P ai 'ticles 
are sufficiently produced. In this case, the increase of Ak due to x-particle 
production is negligible, and resonance stops when q drops down under unity 
(solid curve). In the second case, when |£| is as large as |£| > 10, the growth of 
X fluctuation stops because of the increase of Ak due to x-particle production 
(dotted curve). In this case, resonance ends when the x field reaches the line 
Ak ~ 3q because few instability bands exist for Ak > 3g. 

FIG. 3: The evolution of (x 2 ) as a function of i for £ > 0, m x — in the case of the 

R 2 inflation model ((a) £ = 3, (b) £ = 5, and (c) £ = 10). We find that the 
parametric resonance occurs by the positive ^-coupling, (X 2 ) increases quasi- 
exponentially and reaches its final value (x 2 )f- After then, it decreases by the 
adiabatic expansion, (x 2 )/ takes the maximal value when £ « 5. For larger 
£ (£ > 100), although the growth rate becomes large, the final value (x 2 )f is 
suppressed by r\ term. 

FIG. 4: The final value of (x 2 ) as a function of £ for m x = in the case of the R 2 

inflation model. When £ is positive, (x 2 )f takes its maximal value (x 2 )f — 
1.74 x 10 -4 when £ w 5. For £ < 0, (x 2 )/ takes its maximal value (x 2 )max ~ 
3 x 10~ 4 when £ rs —4 and £ w —40, which is larger than the £ > case. 

FIG. 5: The evolution of (x 2 ) as a function of t for £ > 0,m x — mm the case of the 

R 2 inflation model ((a) £ = 5, (b) £ = 10, and (c) £ = 20). The growth of (x 2 ) 
is suppressed by the mass of x-particles, especially when £ is small. However, 
massive x-P ar ticles of the order m x ~ 10 13 GeV can be created sufficiently 
when £ > 10. 

FIG. 6: The final value of (x 2 ) as a function of m x for £ = 100 in the case of the R 2 

inflation model. (x 2 )f is a quadratically decreasing function of m x . x-particles 
whose mass is m x > 10m ~ 10 14 GeV can not be produced, as is found by 



Eq. (3.12) 



FIG. 7: The evolution of (x 2 ) as a function of i for £ < 0, m x = in the case of the R 2 

inflation model ((a) £ = — 1, (b) £ = —3, and (c) £ = —5). Since the value of \q\ 
is larger as compared with the same absolute value of £, parametric resonance 
occurs even if £ = — 1. The value of (x 2 )/ = 2.63 x 10~ 4 in the £ = —5 case is 
larger than the maximal value (x 2 )f = 1-74 x 10~ 4 in the positive £ case. In 
the negative £ case, (x 2 )/ takes maximal value (x 2 )/ ~ 3 x 10~ 4 at £ « —4 
and £ w -40. 

FIG. 8: The evolution of (x 2 ) as a function of i for £ < 0,m x — m in the case of the 

R 2 inflation model ((a) £ = —5 and (b) £ = —10). When £ = —5, the mass 
effect significantly suppresses the final abundance. However, when £ = —10, 
the mass effect of the order m x = m does not play a relevant role. We find 
that massive x _ P ar ticles of order m x ~ 10 13 GeV can be created sufficiently 
when £ < — 10. 
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FIG. 9: The evolution of (% 2 ) as a function of i for £ > 0, m x — in the case of the 

R 4 inflation model ((a) £ = 10, (b) £ = 30 and (c) £ = 50). We find that 
the parametric resonance occurs by the positive ^-coupling. Since the initial 
fluctuation is small, the final value of (x 2 ) is rather small for £ = 10. However, 
in the case of £ = 30 and £ = 50, the final value of (x 2 ) reaches the value 
(x 2 )/ ~ 10 -5 . For £ > 25, (x 2 )/ does not increase because of the suppression 
effect due to x-particle production. As a result, in the £ > case, (x 2 )/ takes 
almost constant value: (1 ~ 3) x 10~ 5 for 25 < £ < 100. 

FIG. 10: The final value of (x 2 ) as a function of £ for m x = in the case of the R 4 

inflation model, (x 2 )/ takes almost constant value (x 2 )/ ~ (1 ^ 5) x 10 -5 for 
20 < |£| < 100. The maximal value of (x 2 )/ is (x 2 )max ~ 5 x 10" 5 at £ w -35. 

FIG. 11: The evolution of (x 2 ) as a function of t for £ < 0, m x — in the case of the 

R 4 inflation model ((a) £ = -5, (b) £ = -20 (c) £ = -50 ). When £ = -5, 
parametric resonance evidently occurs, but the final value of (x 2 ) is rather 
small. However, in the case of £ < — 20, x-particles are sufficiently produced. 
For -100 < £ < -20, (x 2 )/ takes almost constant value (x 2 )/ « (2 ~ 5) X 10 -5 . 

FIG. 12: The evolution of (x 2 ) as a function of t for £ = — 50,to x = m^ in the case of 

the R 4 inflation model. As compared with the massless case, (x 2 ) is strongly 
suppressed by the mass effect. Generally, massive x - P ar ti c l es which do not 
satisfy the condition Eq. (4.13) can not be created. 
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Fig.3(b) 1=5, m,=0 
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Fig. 4 \ - < % 2 > relation 



v 

a 



-3 

-4 
-5 
-6 

-7 



-100 



-50 







50 



100 



35 



V -8 



-10 - 

-12 ^ 
-14 



Fig.5(a) £= 5, m x =m 



: i , - ; 



5 10 15 20 25 30 

t 



V 



-6 

-8 
-10 
-12 
-14 



Fig.5(b) 1= 10, m ( =m 





- jw.lL.. 


tt!|tl"l'f VvWVw — ^A^WWWW 


: jf 


; 


v» ; 


; 




" 



10 



15 
t 



20 



25 



30 



Fig.5(c) ^=20, m =m 



v 

a 




:-!(» 



Fig. 6 m x - <X > f relation 



N 

V 

a 




10 



12 



rrw 



37 



A 

•V 

s 



-2 

-4 

-6 

-8 

-10 

-12 

-14 
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